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We consider the optimal control problem of transferring population between states of a quantum 
system where the coupling proceeds only via intermediate states that are subject to decay. We pose 
the question whether it is generally possible to carry out this transfer. For a single intermediate 
decaying state, we recover the Stimulated Raman Adiabatic Passage (STIRAP) process which we 
identify as the global optimum in the limit of infinite control time. We also present analytical 
solutions for the case of transfer that has to proceed via two consecutive intermediate decaying 
states. We show that in this case, for finite power the optimal control does not approach perfect 
state transfer even in the infinite time limit. We generalize our findings to characterize the topologies 
of paths that can be achieved by coherent control under the assumption of finite power. If two or 
more consecutive states in an TV-level chain are subject to decay, complete population transfer with 
finite-power controls is not possible. 



PACS numbers: 

I. INTRODUCTION 

Stimulated Raman adiabatic passage (STIRAP) 
achieves coherent population transfer in three-level atoms 
or molecules despite the short lifetime of the intermedi- 
ate level [T]. The key is the creation of a dark state 
produced by overlapping pump and Stokes pulses in a 
counter-intuitive sequence. In the adiabatic limit, the 
intermediate state then never gets populated. STIRAP 
was first demonstrated two decades ago [2] but it contin- 
ues to enjoy great popularity due to its simple, yet robust 
character [3111]. 

Inspired by this, we consider the general A^-level sys- 
tem with a subspace that is free of relaxation, as shown 
in Fig. [1] We consider the case that there is no direct 




FIG. 1: (Color online) A general A'^-level system where the 
states in the subspace S are connected through intermediate 
states outside of 5*. 
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coupling between states in the relaxation-free subspace, 
but the states are coupled by intermediate states that un- 
dergo relaxation. The question we ask here is what kind 
of coupling topology can ensure state-to-state controlla- 
bility on the relaxation-free subspace? We show that this 
question can be reduced to asking what kind of coupling 
topology can ensure unit efficiency of population transfer 
for any two eigenstates in the relaxation-free subspace. 

Our work is closely related to previous studies of STI- 
RAP in multi-level chains [SHE] which showed that un- 
der certain assumptions the dark state condition can be 
generalized from the three-level to the A^-level case. In 
these studies, the decay from intermediate levels was not 
explicitly taken into account. Here, we include the dis- 
sipation. Moreover, we use an analytical formulation of 
optimal control theory which allows us to draw striking 
general conclusions about A^-level systems. In particular, 
we are able to show a relationship between controllabil- 
ity and connectivity of decaying states in the chain: A 
relaxation-free subspace is controllable on the pure state 
space if and only if any two eigenstates in the subspace 
can be connected by a path that never visits two con- 
secutive states that both suffer relaxation. This coupling 
topology includes degenerate levels provided that a gen- 
eralized Morris-Shore transformation exists to replace the 
coupled multi-level system by a set of two- and three-level 
systems and single dark states [16]. 

The problem of controllability in a relaxation-free sub- 
space is closely related to fault tolerant quantum comput- 
ing and decoherence-free subspaces in quantum informa- 
tion science. Given that the Hamiltonian obeys a certain 
symmetry, two or more physical qubits can be employed 
to encode one logical qubit that is free of decoherence [9] . 
The condition for a dark state ensuring STIRAP turns 
out to be equivalent to the condition for a decoherence- 
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free subspace to exist [5] . While in principle it is possible 
to construct quantum gates that preserve the structure 
of the decoherence free subspace [TU] , these gates are gen- 
erally difficult to implement in practice for the following 
reason: the gate operations need to be carried out with 
controls that act on the physical qubits and this intro- 
duces couplings to the decohering subspaces [TT]. This 
raises the question of whether losses can still be avoided if 
the controls are chosen in an optimal way. Previous work 
has discussed whether optimal control can find STIRAP- 
like solutions in A^-level chains [Sl I^UH^ . However, none 
of these studies took the dissipation explicitly into ac- 
count. 

Our paper is organized as follows. In section |llj wc 
study the three-level system. Using optimal control the- 
ory we show that a STIRAP-like process represents the 
globally optimal process for the population transfer in 
the relaxation-free subspace. This STIRAP-like process 
is the infinite time limit of analytic solutions we find for 
the problem with finite time. In section III we study 



a four-level chain system where the two intermediate 
states are subject to decay. We show that with limited 
pulse power, it is not possible to achieve complete pop- 
ulation transfer. Section [IV] generalizes these results to 



A''-level chains. This generalization is used in section [V| 
to state the conditions on state-to-state controllability on 
the relaxation- free subspace. Section [Vl| concludes. 



II. THREE-LEVEL SYSTEM 

We consider a three- level A-system with states |2), 
|3) where |2) suffers relaxation loss with rate k. It is 
well known that population transfer from |1) to |3) is 
possible by STIRAP without populating state |2), i.e., 
complete transfer is achieved in the adiabatic limit [1]. 
In this section, we formulate the population transfer as 
an optimal control problem, taking the dissipation ex- 
plicitly into account. Previous work has addressed this 
problem using a numerical density matrix optimization 
[12j. Here, we formulate the problem analytically using 
the Hamilton-Jacobi-Bellman method which has the ad- 
ditional advantage of allowing us to determine the global 
optimum. We will show that STIRAP arises naturally as 
the solution to the optimal control problem in the adia- 
batic limit. By yielding an upper bound for the transfer 
efficiency in finite time, this formulation also gives some 
insight into the three-level system in the non-adiabatic 
regime. 

The dynamics of the three-level system are described 
by the following effective Schrodinger equation. 



(1) 



where we assume that the detuning is zero. Qp and fig 
are half the Rabi frequencies of the pump and Stokes 
pulses [25], and k is the decay rate of state |2). We want 
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to optimize the transfer of population from |1) to |3) 
within a given time T, i.e., to steer the system from the 
initial state |0(O)) = (1, 0, 0) to the final state such 
that 1 0:3 (T) I is maximized. 

We first make a change of variables, setting xi ~ x'^, 
X2 — ix'2, x^ — ~x'^. The dynamics becomes 



(2) 



Under these dynamics, if we start from the initial state 
(1,0,0), all the state variables will remain real, i.e., this 
change of variables is motivated by the structure of the 
Schrodinger equation and the fact that we are looking for 
state-to-state control. Initial states such as (e'''', 0, 0) can 
be written as e**^ (1,0,0), and since quantum operations 
act linearly on the states, these cases can be reduced to 
the (1, 0, 0) case. The goal is to transfer population from 
(1,0,0) to the final state such that a;3(T) is maximized 
under the controls fip and Q,s[2Q . 

In a realistic setup, both pump and Stokes pulses are 
limited in amplitude, but we will relax this condition: we 
assume that ^Ig is bounded in amplitude by A while fip 
is unbounded. This assumption enables us to solve the 
problem analytically and yields an upper bound on the 
transfer efficiency. Since the amplitude of the pulses is 
usually quite large compared to the relaxation rate, these 
upper bounds are quite tight. We will also show that in 
the adiabatic limit, the condition of unbounded Op can 
be relaxed and STIRAP-like pulses arise naturally from 
the solution of the optimal control problem. 



A. Optimal solution 



To find the optimal pulses, we make another change 
of variables, setting ri — \/ x\ + X2, r2 — X3, and 
tan6' = — , where 9 G [0, The dynamics for (ri,r2) 



are derived from Eq.([2|, 



d 
~dt 



~kcos^ 9 
f2„ cos 9 



-Vis cos 9 




(3) 



Note that fip is now contained in 9 and ri is related to 
the bright-state amplitude in STIRAP. This change of 
variables is not intuitive but crucial for obtaining an an- 
alytical solution. For equations of motion linear in the 
control, the optimal control problem typically becomes 
singular and no conditions to determine an analytical so- 
lution are obtained. After this change of variables, the 
equations of motion are non-linear in one of the controls, 
cos 9. This will allow us to obtain a non- vanishing con- 
dition determining the optimal solution when applying 
Pontryagin's maximum principle. 

Considering the physics of the problem we find that O^ 
should take on the maximal amplitude A throughout the 
process: We start from the initial state (ri,r2) = (1,0) 
and want to maximize r2(T). The population transfer 
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between ri and r2 depends on the rotation speed which 
is determined by flgCosO. At the same time, the popu- 
lation transfer is compromised by ri undergoing decay. 
The effect of decay on ri can be decreased by lowering 
the value of cos9. In order to keep the rotation speed 
between ri and r2 constant, fls needs to be increased. 
For a given rotation speed, the minimum value of cos 6 
and thus the minimum effect of the decay is obtained for 
fls taking its maximum value, fig = A. 

We are now left with determining flp, or, equivalently, 
the angle 9. Defining u = cos 9, we rewrite the dynamics 
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To determine the optimal control we use the prin- 

ciple of dynamic programming |13| and solve for the 
maximum achievable value of r2{T) for all initial points 
(^ij''2)- Starting from (ri,r2), we denote the maximum 
achievable value of r2 by V{ri,r2,t), also called the opti- 
mal return function for the point (ri, r2) at time t. Note 
that for finite time problems, T < cx), the optimal re- 
turn function has an explicit dependence on time, which 
has to satisfy the well known Hamilton-Jacobi-Bellman 
equation, 

dV 

max H{u) = , 



at 



where 



H{u) = ( 



dV 
dri 



dV 
dr2 
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(5) 



(6) 



is the Hamiltonian of the optimal control problem. 

By solving the Hamilton-Jacobi-Bellman equation, we 
obtain the optimal solution to the control problem. The 
detailed derivation is presented in Appendix]^ Here we 
just state the results we need for the subsequent analysis. 
For control time T longer than a critical time Tm the 
optimal control has two distinct phases: 

f ^ for t e [0, t1 

^ ^ \ 1 forte [t,T] ^ ' 

Expressions for Tf^j and t are given in Appendix [X] 
The second stage of this solution, u*{t) — 1, is a result 
of the artificial choice of constraints in our formulation 
(bounded Q.s and unbounded fip). But this need not 
concern us: as we show in Section [IlB| the second phase 
of the solution vanishes in the limit T — >■ oo (adiabatic 
limit), in which case we recover the STIRAP solution. 



B. Recovery of STIRAP 

We now show that in the limit T ^ oo, our optimal 
pulse corresponds to STIRAP. The optimal Rabi frequen- 
cies are derived from ri, r2 and u* . From Eq.([2]), it fol- 
lows that 



dt 



X2 = ripXi — kX2 — ^sXz ■ 



(8) 
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FIG. 2: (Color online) Optimal half Rabi frequencies Q,s{t) 
and Q.p{t) for finite control times T and relaxation rates k 
(with the amplitude bound on Q,s set to one, A = 1). 



Substituting X2 
obtain 



riu, xi = ri 



and X3 — r2, we 



rip 



—ku^ri — flsUr2 + riii + kriu + ^sT2 

T\\J\ — u^ 



(9) 



When the control time goes to infinity, T — > oo, the 
switching time, r, also approaches cxd. This follows from 
T—T being smaller than the critical time Tm as explained 
in Appendix [Xj If r becomes very large, then according 
to Eq. Q, w*(i) is very small. Therefore X2 = r\u* ^ 0, 
i.e., the level |2) is not populated. Substituting w ~ 
and ii ^ into Eq. (|9| , we see that 



n 



X3 



= '-In. 



(10) 



So at each time point, ^2 ~ 0, ^ = f^, i.e., the sys- 
tem is in the state a;i|l) + xajS), where = which 
corresponds to the dark state in STIRAP. 



Note that in this limit, an upper bound can also be 
assumed for Qp in addition to the one for ilj. Since it is 
the ratio Qp/Q. that matters, one can simply lower ils if 
fip exceeds its upper bound to maintain the ratio Op/fis. 
For example, one could assume the same bound A on 
flp. In the optimal solution shown in Fig. [21 one would 
then rescale time whenever flp hits A. This obviously 
avoids the shoot-off of ftp to infinity visible in Fig. [2] at 
late times. Such a rescaling of amplitude corresponds to 
changing the unit of time: When Qp > A, a, new time 
variable dr — 



-dt is defined. In this new unit of time. 
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the dynamics become 




(11) 



The relaxation rate k does not affect the dynamics in the 
infinite time hmit, so k and ^k are effectively the same, 
and the ratio of the Rabi frequencies still satisfies the 
optimal condition in Eq. (10). Since the reparametrized 
time remains infinite in the limit T — >■ oo, optimality of 
our solution is not affected by changing the unit of time. 

One might wonder in this case where the characteristic 
time delay between and fip is hidden in our solution. 
The point is that STIRAP is determined by the over- 
lap of the pulses. The rising part of the Stokes pulse 
when the pump is zero and the falling part of the pump 
pulse when the Stokes is zero do not affect the system 
dynamics. Our solution only contains the crucial over- 
lapping part, similar to the "shark-fin" pulses discussed 
in Ref. [H]. The Stokes pulse starts out fairly flat at 
the upper bound and falls down as the pump pulse is 
rising up to the bound. So the delay between the pulses 
in standard STIRAP corresponds to the time the pump 
pulse takes to rise to the upper bound. In our solution, 
a rising edge of the Stokes pulse and a falling edge of 
the pump pulse could be added if one wishes to obtain a 
more realistic pulse shape. 

To summarize the similarities and differences with the 
conventional STIRAP solution, we drop the constraint on 
the bound of the pulse amplitude fJp and obtain an ana- 
lytical solution for the optimal control (cf. Eqs. ([7|,([9])) 
in finite time. While in general the optimal pulse shape 
can be found only numerically, for infinite time we ob- 
tain a completely analytical solution (cf. Eq. (lOl). By 



solving the Hamilton- Jacobi-Bellman equation, we have 
proven that this solution is the global optimum. The rise 
of the Stokes pulse and the fall of the pump pulse are 
missing. These portions of the pulses are irrelevant in the 
infinite time limit. Thus, our analytical solution confirms 
that the essential feature of STIRAP is the time-ordering 
of the pump and Stokes pulse where they overlap. 



FIG. 3: (Color online) A four-level system: Population shall 
be transferred from |1) to |4) via the intermediate states |2) 
and 1 3) which suffer relaxation. 



complete transfer with limited power even if we wait in- 
finite time. 

The dynamics of this system are described by the fol- 
lowing effective Schrodinger equation. 
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where ftp and ilg are the Rabi frequencies of pump and 
Stokes pulses, 17/ is the Rabi frequency of a pulse cou- 
pling 1 2) and |3), and k denotes the decay rate. As in the 
previous section, we assume zero detunings. Since we 
are only interested to find out whether there are schemes 
which avoid populating the intermediate states, the exact 
value of the relaxation rate is not important here. For 
simplicity, we assume that the intermediate states suffer 
the same amount of relaxation. It is straightforward to 
generalize to the cases with different decay rates. 

We want to find the optimal way to transfer population 
from |1) to |4) within a given time T, i.e., the optimal 
way of steering the system from the initial state \<p{0)) = 
(1,0,0,0) to the final state \(j){T)) such that \x'^{T)\ is 
maximized. 

Analogously to Section In) we make a first change of 



variables, letting xi 



-i5 X2 — 1X2, x^ 
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The dynamics become 
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(13) 



III. FOUR-LEVEL SYSTEM 

In this section, we extend our method to a four-level 
chain system as shown in Fig. |3] Again we explicitly in- 
clude the decay from the intermediate levels. We can thus 
demonstrate that the four-level chain differs fundamen- 
tally from the three-level system where STIRAP-like pro- 
cesses can transfer the population fully in the adiabatic 
limit. It turns out that in four-level systems with two 
decaying intermediate states it is not possible to achieve 



The state variables are now all real numbers if we start 
from the initial state (1,0,0,0). We want to transfer 
from (1, 0, 0, 0) to the final state such that Xi{T) is max- 
imized under the controls fip, VLi and fi^ at given time 
T. We again relax the control constraint, assuming that 
VLi is bounded in amplitude by A, but fJp and ^Is are 
not bounded. We show that even with these relaxed con- 
trol constraints, it is not possible to achieve unit transfer 
efficiency. 

To solve the problem, we make a second change of vari- 
ables, letting ri = \/x\ + x\, r2 = \/x\ + x% tan^i = 
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d_ 
~dt 



tan 9i 

r2 



X4 
X3 ' 



The dynamics of (ri,r2) become 



— fccos (?i — 0/ cos 6'i cos ^2 1 ^1 
17/ cos 01 cos 6*2 — fccos^6'2 / V ''2 

(14) 

This looks familiar since we have almost the same equa- 
tion as for three- level system, cf. Eq.([3]). We first ob- 
serve that flj should always take on the maximal ampli- 
tude A: if rij < A, we can always increase it to A while 
lowering cos^i and cos 6*2 such that the rotation speed 
ri/ cos 9i cos 6*2 remains the same, but the effect of decay 
on ri, r2 is decreased. The problem therefore reduces to 



d_ 
~dt 



—ku\ 

AuiU2 
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(15) 



where ui = cos 9i and U2 — cos 62 ■ The same dynamics 
arise in Nuclear Magnetic Resonance and analytical so- 
lutions to this control problem were obtained using the 
optimal control technique in Ref. [14]. We describe the 
characteristics of the optimal pulse sequences and refer 
to Ref. [13] for more details. 



Case I: If T < 



where C = then u\(t) 



^2(0 — 1 throughout. That is, we obtain a hard pump 
pulse at t = 0, flipping the angle 9i by it/2, i.e. trans- 
ferring all population from xi to X2, and a hard Stokes 
pulse at i = T, flipping 62 by 7r/2, transferring popu- 
lation from to x^. At intermediate times, the op- 
timal Rabi frequencies for pump and Stokes pulses are 
zero. The efficiency of the population transfer, rj^, is ob- 
tained by integrating Eq. ( 15 ) with the optimal controls 
ul{t) = ul{t) = 1, 



riT = exp[-fcr]sin(v4T) . 



(16) 



This is smaller than unity for all T > 0. 

Case II: For larger control times, T > SSL^Ml^ ^ jg 
not optimal to put all population immediately into the 
decaying level X2- The optimal trajectory then has three 
distinct phases. For < t < r, where r is a function 
of T, U2 {t) = 1 and ul (t) is increased gradually from a 
value uj(0) < 1 to ^^(t) = 1. This corresponds to flp{t) 
rising from its initial value r2p(0) to infinity while fls{t) 
remains zero throughout the first phase, cf. Fig. [4] In 
the second phase, for time t < t < T — t, the optimal 
controls are ul{t) = U2{t) — 1. This corresponds to 
both r2p(i) and r2s(0 being zero. Finally, for t > T — r, 
ul{t) = 1 and m5(*) is decreased from U2(T — r) — 1 
to U2{T) = Mi(0). fls{t) is thus decreased from infinity 
to its final value ris(T) which is equal to r2p(0), while 
^p{t) remains zero. The parameter r determining the 
switching times is calculated from the following equation. 
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FIG. 4: (Color online) Optimal half Rabi frequencies Qp{t), 
fli{t), and 0,s{t) (Case 2) for finite control times T with re- 
laxation rate fc = 1 and the amplitude bound on Qi set to 
one, A — 1. For very short T (Case 1), the optimal Rabi fre- 
quencies f2s(t) and Qp{t) correspond to instantaneous pulses 
of infinite power at t = and t = T; for T — >■ 00, the hold pe- 
riod with zero ^ls{t), rip{t) in the middle of the time interval 
disappears. 



and 



k(t) = 1 + 2^^ 



2CVl + coth [A^/l + S,^ r + 2sinh"^C 



In the limit of infinite time, r = -j; in this case the 'hold' 
phase in the middle, where only il/(t) is non-zero, dis- 
appears (cf. the blue and green curves in Fig. |4]) . The 
optimal solution thus corresponds to the intuitive pulse 
sequence of pump first, then Stokes, not a STIRAP- 
like solution characterized by a counter-intuitive pulse 
sequence. 

The efficiency of the population transfer, rjT, in case 2 
is expressed for finite time T in terms of the angles as 



exp(g(7i -72))(1 -esin2g2) 
sin (71 + 72) 



(18) 



In the limit that T goes to infinity, 71 = 72 = 
tan^^ ^-^/l + — , and rjx approaches i], the maxi- 
mum transfer efficiency given by 



V 



(19) 



We see that the transfer efficiency can reach unity only 
for ^ = 0, i.e., ^ = 0. This is the case of infinite power, 
i.e., the decay is much smaller than the maximum Rabi 
frequency coupling the intermediate states. When only 
limited power is available, the transfer efficiency is always 
less than unity. 

Even for the case of a multiphoton resonance, no an- 
alytical optimal solutions are known for non-zero detun- 
ings. As shown in Refs. [3 [5], adiabatically eliminating 
one of the two decaying intermediate levels allows one to 
recover a three-level system and thus the STIRAP solu- 
tion. However, perfect adiabatic elimination in the pres- 
ence of decay requires infinite power of the field coupling 
the intermediate levels. We believe that our conclusion 
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for a finite power transfer efficiency of less than unity 
holds also for non-zero detunings. 

To support this claim, consider the four-level system 
with non-zero detunings A2, A3, 
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Assume complete population transfer is possible in this 
system. Clearly, for this to hold the population of the 
intermediate states must remain zero during the process. 
If this is the case, the values of the detunings would not 
affect the transfer and we can therefore replace the de- 
tunings by zero. However, in this case the system re- 
duces to Eq. ( 20 1 , which we showed to be uncontrollable 
leading to a contradiction with our 
Thus we claim that even 



with finite power, 
assumption of controllability, 
with non-zero detunings, complete population transfer is 
not possible with finite pulse power if there are multiple 
decaying intermediate states. Since detuning relative to 
the intermediate states is known not to affect STIRAP 
efficiency [T] , this result is fully consistent with our prior 
expectations. 



IV. GENERALIZATION TO A-LEVEL SYSTEMS 

In the previous section we showed that with two con- 
secutive intermediate decaying states, the transfer efh- 
ciency with finite power is always less than one. If more 
than two intermediate decaying states are present, the 
transfer efficiency gets even worse. For example, suppose 
we have a chain of five states with three intermediate 
decaying states. 
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With a change of variables, letting xi = x[, X2 



X3 
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Xi ~ —ix'^ and x^ — Xg, the dynamics become 
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Introducing new variables y^, the fifth state and fourth 
state can be combined, 
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This reduces the system to an effective four-state system, 
/ 

V 

(21) 

where tamO — |j. The transfer efhciency for the dy- 
namics of t/i, . . . ,?/4, represents an upper bound to the 
transfer efficiency for the dynamics of . . . ,^5, since 
2/4 > 2:5 : If a control scheme transfers an amount of pop- 
ulation, 77, from xi to then this scheme also transfers 
at least an amount of population 77 from yi to yi. Next 
we show that the efficiency of transfer from yi to y^ ac- 
cording to the dynamics of Eq. (21 ) is upper bounded by 



the transfer efficiency for the dynamics of Eq. (131 of the 
four-level chain. To make the comparison transparent, 
we consider the following dynamics. 
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(22) 

Clearly the efficiency of transfer from yi to j/4 is upper 
bounded by the efficiency of transfer from y'l to 2/4 since 
the only difference between these two dynamics is that 
2/4 is subject to decay while 2/4 is not. The efficiency 
of transfer from y'l to 2/4 in turn is upper bounded by 
the transfer efficiency for the four-level chain, Eq. (13 1, 



since the controls in Eq. (22 1 are more restricted than 
the controls in Eq. (13). This can be seen as follows: If 

reaches an 
(|T3l, 



(221 



2/1 = a;i , 2/2 = a;2 , 2/3 = 2:3 , 2/4 



a control scheme for the dynamics, Eq. 
efficiency 77, then simply setting fl^ — ^2 cos d in Eq. 
the same transfer efficiency ry is obtained for the four-level 
chain. The inverse step of setting ^2 = is, however, 
not always possible, since this may lead to infinite ^2- 

From this line of argument, we see that the transfer 
efficiency of the five-state chain with three intermediate 
decaying states is upper bounded by the transfer effi- 
ciency of the four-state chain with two intermediate de- 
caying states. Analogously, we can show that the transfer 
efficiency of the A^-state chain ( A^ > 4) with N —2 inter- 
mediate decaying states is upper bounded by the transfer 
efficiency of the A^-state chain with A' — 3 intermediate 
decaying states. That is for chains where all intermedi- 
ate states are subject to decay, the transfer efficiency is 
monotonically decreasing with increasing length of the 
chain. In particular, we have the interesting result that 
for any chain with two or more consecutive intermediate 
decaying states, the efficiency will be less than unity. 

Our analytical results allow us to draw a general in- 
ference about previous work on population transfer in 
A^-level systems. Extensions of STIRAP from the three- 
level system to multi-level chains have been investigated 
since the early days of STIRAP (see Ref. [1] and refer- 
ences therein). Although these mechanisms are designed 
to keep the population in the intermediate levels as small 
as possible, close inspection reveals that none of these 
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succeed in completely avoiding intermediate population 
for finite power pulses, even for T — > c». In particular, 
we note that none of these schemes can keep the pop- 
ulation for two consecutive intermediate states at zero. 
We illustrate this with three of the generalized STIRAP 
schemes. In Ref. [5], a generalization of STIRAP for A^- 
level chains was proposed for N odd. In this scheme, 
the population in all the intermediate even levels can be 
kept at zero but placing population in the intermediate 
odd levels cannot be avoided. Clearly it is impossible 
to keep two neighbouring states empty. In Ref. [6 , a 
STIRAP-like solution for A^-level chains was found nu- 
merically, and termed straddling STIRAP. It consists in 
choosing the Rabi frequencies coupling the intermediate 
states to be at least one order of magnitude larger than 
fls and rip and overlapping in time with both fls and 
flp. Inspection of the solution reveals that all interme- 
diate levels acquire some population, and therefore it is 
impossible to keep two consecutive levels empty without 
using infinite power. The straddling STIRAP was an- 
alyzed further both analytically and numerically [71 [8] . 
It was clarified that for N even, a non-zero detuning of 
the lasers coupling the intermediate states is required to 
obtain a STIRAP-like solution while for zero detunings, 
an intuitive (Rabi) pulse sequence is found. Moreover, 
it was shown that in the dressed state picture, a very 
strong coupling between the intermediate states moves 
the intermediate states out of resonance such that they 
are decoupled and effectively a two-level system {N even 
with zero detuning) or a three-level system (N odd or N 
even with non-zero detuning) are recovered. It can be 
seen that whether TV is odd or even, avoiding population 
in two consecutive levels requires infinite power. Taking 
dissipation explicitly into account, it is clear that if two 
consecutive intermediate levels are subject to decay, unit 
transfer efficiency is impossible at finite power. Our an- 
alytical results including dissipation are consistent with 
this analysis. 



V. CONTROLLABILITY IN THE 
RELAXATION-FREE SUBSPACE 

In this section, we will generalize the results of 
the previous sections to state-to-state controllability on 
relaxation-free subspaces S. Based on the results of the 
previous sections, we will characterize the relaxation-free 
subspaces that are finite-power controllable on the pure 
state space. Our argument is based on the assumption 
of selective control, i.e. the assumption that Qs, and 
any other coupling, can be tuned independently. This 
corresponds to the bare, field-free Hamiltonian having no 
degenerate levels or frequencies. With selective control, 
controllability on the relaxation-free subspace becomes 
equivalent to connectivity |15j . 

Our central result is that a relaxation-free subspace is 
controllable on the pure state space if and only if any two 
eigenstates in the subspace can be connected by a path 




1 2 

b 

a 



FIG. 5: (Color online) Examples for (a) a system that is 
controllable in the subspace {|1), |2), |3), |4)} since any pair 
of eigenstates can be connected by a path that never visits 
two consecutive states that suffer relaxation; (b) a system 
that is not controllable in the subspace |2), |3), |4)} since 
one has to pass through three consecutive states that suffer 
relaxation to connect the states |1) and |4). 




1 2 



FIG. 6: (Color online) The eigenstates |1) and |4) are con- 
nected by concatenating paths of type (I) and type (II) . 



that never visits two consecutive states that both suffer 
relaxation. This is illustrated in Fig. [5] with part (a) 
displaying an example where any pair of eigenstates in 
the relaxation-free subspace 5* = {|1), |2), |3), |4)} can 
be connected by a path that never visits two consecu- 
tive states suffering relaxation, i.e. this system is state- 
to-state controllable in the space S = {|1), |2), |3), |4)}. 
The system shown in Fig. [5]jb) is not controllable in the 
relaxation-free subspace S = {|1), |2), |3), |4)} since one 
has to pass through three consecutive states that suffer 
relaxation in order to connect the states |1) and |4). 

Clearly the condition is necessary. As from the previ- 
ous section, we know that if any path has two intermedi- 
ate states outside of the relaxation-free space, complete 
population transfer is not possible. Hence the system is 
not controllable. The condition is also sufficient. If there 
are never two states in a row that suffer relaxation, the 
control found in Section |ll] allows us to traverse one in- 
termediate relaxing state without losses. Concatenating 
such processes gives the result. 

We can connect the controllability condition with the 
coupling topology of the system: The condition is ful- 
filled if and only if the eigenstates in the relaxation-free 
subspace are (I) directly connected via paths in the sub- 
space; (II) connected by one intermediate state which 
suffers relaxation; (III) connected by concatenations of 
paths of type I and type II as sketched in Fig. [6j Note 
that this coupling topology includes degenerate levels in 
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the relaxation-free and the relaxing subspaces provided 
that a generalized Morris-Shore transformation can be 
employed to replace the coupled multi-level system by a 
set of two- and three-level systems and dark states (single 
levels) [H]. 

We first show that with this coupling topology we can 
achieve coherent population transfer between any two 
eigenstates. Controllability is then achieved by applying 
sequences of these operations. First, if two eigenstates 
are connected via (I) we can obviously achieve coher- 
ent population transfer between them. It remains to be 
shown that coherent transformation is possible for two 
eigenstates that are connected by one intermediate state 
suffering relaxation. This is achieved by combining our 
results of Section [IT] with the fractional STIRAP devel- 
oped by Vitanov et al. [ 17) to generate arbitrary coherent 
superpositions of |1) and |3) from the initial state 
For two eigenstates |1) and |3) that are connected by 
an intermediate state |2) suffering relaxation, a coherent 
transformation from state |1) to cos/3|l) — e*'^sin/3|3) is 
implemented by (i) adding a phase </> to the Stokes pulse 
such that the equations of motion read 

^ /x[\ / \ /x[\ 

- a;U = U^p -ik flse^'^ h^2 ' (23) 
V ^3 / V n^e^^ J \x'J 

and (ii) varying ^ adiabatically from to tan/3. A 
general coherent transformation can be implemented by 
first applying the time-reversed version of this procedure 
in order to transfer the initial state to |1) and then to 
transfer |1) to the desired final state with the control 
scheme of Section [llj We have thus shown controllability 
for eigenstates connected by (II). Obviously the results 
for (I) and (II) can be combined which yields controlla- 
bility for (III) completing the proof of sufficiency. 

So far we have considered only zero detunings. How- 
ever, the result of controllability on the relaxation-free 
subspace also holds for non-zero detunings that can be 
represented by complex values for the decays. To see 
this, consider first a three- level system. It is well-known 
that the values of the detuning and the relaxation rate 
of the intermediate state do not affect the STIRAP effi- 
ciency [T]. This argument carries over to general N-level 
systems with non-zero detunings of the decaying states. 
If the system is controllable it can be viewed as a con- 
catenation of two or more three-level STIRAP systems 
(and possibly type (I) paths), and then the detuning is 
irrelevant. If the system is uncontrollable, the detuning 
cannot make it controllable, based on the argument pre- 
sented at the end of Section Hill 

VI. DISCUSSION AND CONCLUSIONS 

We have considered the optimal control problem of 
transferring population in a quantum system between 
states in a subspace free of dissipation, where the trans- 
fer has to proceed via states that are subject to decay. 



We treated only the case of resonant controls with fixed 
carrier frequency and phases, controlling only the ampli- 
tudes of our pulses as a function of time. Such situations 
occur frequently in atomic and molecular physics appli- 
cations. For example, transfer between different levels in 
the electronic ground manifold can proceed via Raman 
transitions. In quantum information applications, stable 
qubit states are often connected via auxiliary states that 
are subject to decay. In particular, this may be the case 
for logical qubits encoded in a decoherence-free subspace. 

We have obtained analytical solutions to this optimal 
control problem by solving the Hamilton-Jacobi-Bellman 
equation for the optimal return function. For a single 
intermediate decaying state, we have recovered the Stim- 
ulated Raman Adiabatic Passage process [TJ as the glob- 
ally optimal solution in the limit of infinite time. Perfect 
state transfer is achieved only in this limit. This is in 
accordance with experimental realizations where at best 
99.5% state transfer were achieved [18 . 

In Ref. [24 , the STIRAP solution in a three-level 
system was previously obtained using geometric control 
methods. There, the optimal return function was specif- 
ically designed to avoid the hard pulses obtained by us 
(which cause the sudden population transfer from level 
1 to 2 at i = 0). Note that generalizing the results of 
Ref. [21] from the three-level system to TV-level systems 
is hampered by the system's state being represented in 
terms of six real variables, compared to two variables, ri 
and 7^2, in our case. 

In contrast to the analytical solutions presented here 
and in Ref. Refs. employ numerical opti- 

mization procedures based on the calculus of variations. 
Our current work may help to clarify the disagreement 
in the literature on whether STIRAP is obtained as a 
solution to an optimal control problem [5J [201 [12] '■ the 
assertion of Ref. 20 that adiabatic passage population 
transfer cannot be obtained as the solution to an opti- 
mal control problem implicitly assumes finite pulse flu- 
ence and finite control time. Yet adiabaticity, strictly 
speaking, does not comply with these assumptions. 

We have also presented analytical solutions for the case 
of population transfer that proceeds via two consecu- 
tive intermediate decaying states. In particular, we have 
shown that the optimal control does not yield perfect 
state transfer even in the limit of infinite time, unless the 
pulse coupling the intermediate levels has infinite power. 
This gives an analytical framework for understanding an 
earlier control solution termed straddling STIRAP that 
was obtained numerically [6 and that is essentially based 
on adiabatically eliminating the intermediate levels ^TjiS]. 
Taking dissipation explicitly into account, we have clari- 
fied that the adiabatic elimination of the decaying levels 
is possible only in the limit of infinite power. 

Finally, we have generalized these results to charac- 
terize the topologies of paths that can be achieved in A^- 
level systems by finite-power controls and in the presence 
of dissipation. Population transfer with unit efficiency 
is only possible if each decaying state is connected to 
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two non-decaying states. Complete population transfer 
is then achieved in the adiabatic limit, i.e., in a sequence 
of STIRAP processes. Finite-power state-to-state con- 
trollability on the relaxation-free subspace is thus equiv- 
alent to connectivity 15j, augmented by the condition 
that only one out of two consecutive levels may be sub- 
ject to dissipation. 
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Appendix A: Optimal control for the three-level 
system 

In this Appendix, we derive the solution of the optimal 
control problem for the three level system. As introduced 
in Section |TT] we are going to solve the Hamilton Jacobi 
Bellman equation, Eq. ([5 ) with the classical Hamiltonian 
of the control problem defined in Eq. Q. Introducing 
adjoint variables Ai, A2, 



Ai = 



dV 

dri 



A2 = 



dV 



(AI) 



where V denotes the optimal return function, the Hamil- 
tonian of the optimal control problem can be expressed 

as 



H{u) 



-kXiriu^ — AXir2U + AX2riu 

-fcAlTi ■ 



A 
2k 



A2 
Ai 



r2 
ri 



-fcAlTl 

A 
2k 



A2 
Ai 



A 
ri 



A2 
Ai 



(A2) 



We rewrite the ratios appearing in Eq. ( A2 ) in terms of 
variables a and b. 



A2 
Ai 



r2 



(A3) 



The optimal return function is a non-decreasing function 
of ri, r2 [starting from a larger ri(0) or ?'2(0), one can 
achieve a larger 7-2 (T)]. Due to Eq. ( |Al[ ) we therefore 
find Ai > 0, A2 > and hence a > 0. Since Ai > and 
ri > 0, maximizing H with respect to u is equivalent to 
minimizing the function 



If a — 6 < 0, then the solution is the trivial one, u* = 0. 
We therefore conclude that obtaining a non-zero control 
requires a — b > 0. Later we will show explicitly that 
a — 5 is a non-decreasing function of time. Since 6(0) — 
and a(G) > 0, we have a{t) — b{t) > for all times t. We 
further distinguish two cases. 

Case I: If ^(a — &) > 1, then the minimum of f{u) is 



achieved at the maximum value that u 
u* = 1. 

Case II: If < 



cos 9 can take, 



f{u) is achieved at u 



(a — fe) < 1, then the minimum of 
* = i:{a-b). 



2k 



It is a standard result that, along the optimal trajec- 
tory (ri(<), r2(t)), the adjoint variables (Ai(t), A2(t)) sat- 
isfy the equations 



dAi 
dt 



dH 
dri 
dH 
dr2 



I.e., 



~dt 



Ai 
A2 



kv? 
Au 



-Au 




Ai 
A2 



(A4) 



with the term inal condition Ai(T) ~ 0, A2(r) ~ 1. 

With Eqs. (|A4| and Q, we can derive the dynamics 
for a and b along the optimal trajectory. 



— a — Aua^ 
dt 



dt 



Au9 



ku^a - 
- ku^b 



Au, 
Au . 



Therefore 



dt 



(a-b) = {a + b)u[A{a - 6) - ku] 



(A5a) 
(A5b) 

(A6) 



For the optimal trajectory starting at (ri,r2) = (1,0), 
b(0) = 0. Depending on a(0) we have the following cases. 
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Case A: If a (0) > ^, we start in Case I, u* 



1. From running it backwards, we obtain 



equation (A6), it follows that 
d 



dt 



{a-b) = {a + b)[A{a - b) - k] > k{a + b) > , 



i.e., a — b is non-decreasing. Therefore we will remain in 
Case I for the whole time interval. Substituting u* — 1 
into the dynamical equation for (Ai,A2), Eq. (A4), and 



1 / Ai(T-<) 
dt [HT-t) 



A 

-A 



with the initial condition Ai(T) 
grating this equation yields 



Ai(r-t) 



(A7) 



0, A2(T) = 1. Inte- 



AiW 
A2(t) 



1 

2u> 



(A8a) 



where 



and 



if fc2 > 

if fc2 < 4^2 



Al(<) 
A2(t) 



3-^('^-*)[A(T-<) + l] 



(A8b) 



for fc2 



4^2 



Solving the equation a(0) = ^ = f 
for T, we obtain a critical time, T/\/. When T < Tm, 
the optimal control is to set u* = 1 for the whole time 
[0,T]. An analytical expression can be derived for Tm- 
For example, when k^ > 4^12^ 



Tm — 



log 



2(A2+2fc2) 
Vfc2 - 4^12 



(A9) 



The other cases can also be worked out. The critical time 
Tm as a function of the decay rate k and the amplitude 
bound A is displayed in Fig. [Tj Tm takes large values 
corresponding to the trivial optimal solution u* = 1 for 
all is [0,T] only for small decay rates and small ampli- 
tude bounds. As k and A increase, Tm quickly becomes 
fairly small (for example, Tm = 0.06 ioT: k = A — 10) and 
the optimal solution is determined according to Eq. ([7]). 

Case B: If a( 0) < ^, we start in Case II, u* — ^{a-^)- 
From Eq. (A6), we get in this case 



dt 



(a-b) 



ia + b)^ia-b) 



A{a-b) (a -6) 



A^ 
4fc 



(a + 6)(a-&)2 >0 



(AlO) 



If a(0) — 6(0) = 0, then it will remain zero for the whole 
time. We will see that this only occurs as T — >■ oo. If 
a(0) > 0, then a{t) — b{t) is strictly increasing, and at 
some time r, it will reach We then switch to Case I, 
setting u* = I afterwards. So in this case, the optimal 
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FIG. 7: (Color online) The critical time Tm as a function of 
decay rate k and amplitude bound A. For control times T 
shorter than Tm the optimal control is simply set to w* = 1 
for all t G [0,T], for control times larger than Tm, Eq. Q 
applies. 



control is u*{t) = ^[a{t) — b{t)] for t < r and u*{t) — 1 
for T < i < T. 

We now show how to c alculate the switching time t of 
Case B. Using Eqs. (A5) and u* — ^{a — b), one can 
show that within the time interval [0, r], 



d 
dt 



' k 



Together with the initial condition, °|o)^fc(o| 
in mind that 6(0) = 0, this yields 



a(r) - 6(t) _ 1 
a(r)+6(r) - 



(All) 



1, keeping 



(A12) 
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Since at time t = t, air) — b{T) = 

a(T) = At + 
b{T) = At. 



= we obtain 
2k 



(Al3a) 
(Al3b) 



In the time interval [t, T], u* = 1, and we can again 



run Eq. (A7) for (Ai, A2) backwards in time from T to r. 



This yield s another expression of a(T) from Eqs. ( A8a 



and jA8b). Setting it equal to At + we obtain the 
cning 



switching time t. 

Note that T-t <Tm- From Eqs. ( |A13[ ), the fact that 
b{t) is non-negative for all times t and a{t) — b{t) is an 
increasing function, it follows that 



a{t) > 



2k 
~A 



for t e [t, T] . 



The assumption T — t > Tm then leads to a contradic- 
tion: If T - r > Tm, then a{T - Tm) = 2k /A at time 
t = T-TM. 

Next we evaluate the value of the optimal control 
u*{t) = ^[a(0 - b{t)] for t < t. We know that in 
the interval [0, t], a~b satisfies the dynamical equation, 
Eq. ( |A10 1, which can be rewritten 



dt 



(a-b) 



A^a + b 
Ak a-b 



{a-bf 



(AM) 



From Eq. (|A12l), we get ^ 



— + 1. Substituting 
it into the above differential equation and solving it, we 
obtain a{t) - b(t) for t G [0,t], 



a{t) - b{t) 



2k 



A 



f^) + 2k{T - t) + I 

So in Case B, the optimal control is obtained to be 

for t e [0, r 



(A15) 



u*{t) 



^A2(r2-t2)+2fe(r-t) + l 



1 



for t £ [r, T] 



(A16) 



In summary, there is a critical time, Tm, which depends 
on the relaxation rate k and the amplitude bound A and 
determines whether the control is switched or not: For 
T < Tm the optimal control is just set to one, u*{t) = 1, 
for all times t e [0,T], and for T > Tm, 

f ^ for t e [0, t1 

[1 for i e [r, T] 

The optimal control u* (t) is thus determined by the sys- 
tem parameters A and k and the switching time r which 
is obtained by matching the dynamics of a{t) at t — t, 
cf. Eqs. (lAlSal) and (|A8|. 



We derive the optimal Rabi frequency, ^p(t), from 
riit), r2{t) and u*{t). From Eq. it follows that 



dt 



X2 



fipXi — kX2 — ^sXz ■ 



(A17) 
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FIG. 8: (Color online) The optimal transfer efficiency, r2{T) 
as a function of the control time T for different relaxation 
rates k (the maximum Rabi frequency, or pulse power, re- 
spectively, is limited by ^ = 1). The curves are fairly con- 
cave, reaching a high transfer efficiency in relatively short 
time. This is compatible with the fact that in real experi- 
ments stimulated Raman adiabatic passage is actually done 
in relatively short time. To reach 100% transfer efficiency, 
however, infinite time is needed. 



Substituting in Eq. ( A17) X2 — riu, xi — riVl — 
r2, we obtain 



and 



X3 



nr. 



-kv^ri — Aur2 + riii + kriu + Ar2 



ri 



(A18) 



It is difficult to obtain a closed form of flp{t), but Fig. [2] 
shows a few numerical examples. The solution u*{t) = 1 
occurs toward the end of the interval [0,r]. For finite 
control times T, this solution for u*{t) corresponds to 
n,p(t) being infinity. For T — >■ 00, a rescaling of time 
leads to finite Op(t) as explained in Section 



II B 



The interpretation of the optimal pulses is as follows: 
For small control time T, the major limitation for the 
population transfer is not due to relaxation, but the lim- 
ited available time. The optimal choice u* = 1 maximizes 
the transfer speed, but also maximizes the decay of ri (t) , 
respectively X2(i), as can be seen from Eq. Q. However, 
for a small available control time T, the gain obtained 
by maximizing the desired transfer at each moment is 
more important than the relaxation losses. As the con- 
trol time T increases, the relaxation degrades the perfor- 
mance more and more and the choice u = 1 ceases to be 
optimal. For finite time the optimal solution becomes a 
compromise between maximizing the transfer speed and 
minimizing the decay. When time goes to infinity, mini- 
mizing the relaxation loss becomes more important than 
the transfer speed. 

Substituting the optimal control u* (t) into the dynam- 
ics, Eq.Q, and integrating it yields the value of r2{T). 
For finite T this gives an upper bound for the maximum 
achievable population transfer due to the possibility of 
r2p(t) going to infinity. As shown in Fig. [s] the upper 
bound approaches unity in the limit T — )• 00 even for 
large decay rates fc, illustrating the recovery of STIRAP. 
For small and moderate decay rates, the convergence to- 
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ward unit efficiency is much faster, reflecting the easier control problem. 
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